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ABSTRACT 



Aims. Following the advent of increasingly sensitive X-ray observatories, deep observations of early-type stars became possible. 
However, the results for only a few objects or clusters have until now been reported and there has been no large survey comparable to 
that based upon the ROSAT All-Sky Survey (RASS). 

Methods. A limited survey of X-ray sources, consisting of all public XMM-Newton observations (2XMMi) and slew survey data 
(XMMSLl), is now available. The X-ray counterparts to hot, massive stars have been searched for in these catalogs. 
Results. About 300 OB stars were detected with XMM-Newton. Half of them were bright enough for a spectral analysis to be possible, 
and we make available the detailed spectral properties that were derived. The X-ray spectra of O stars are represented well by low 
(<lkeV) temperature components and seem to indicate that an absorption column is present in addition to the interstellar contribution. 
The X-ray fluxes are well correlated with the bolometric fluxes, with a scatter comparable to that of the RASS studies and thus larger 
than found previously with XMM-Newton for some individual clusters. These results contrast with those of B stars that exhibit a 
large scatter in the Lx - Lbol relation, no additional absorption being found, and the fits indicate a plasma at higher temperatures. 
Variability (either within one exposure or between multiple exposures) was also investigated whenever possible: short-term variations 
are far more rare than long-term ones (the former affects a few percent of the sample, while the latter concerns between one third and 
two thirds of the sources). 

Conclusions. This paper presents the results of the first high-sensitivity investigation of the overall high-energy properties of a sizable 
sample of hot stars. 

Key words. X-rays: stars - Stars: early-type 



1. Introduction 

Soon after the discovery of X-ray emission from hot stars 30 
years ago, it was proposed that a correlat ion exists betw e en the 
X-ray and bolometric luminosities (e.g. iHarnden et al.L 119791: 
|PaUaviciniet"d].ll98lh . of the form L™^''' ~ 10"^ x Lboi- This 
relation probably reflects the dependence of the X-ray emission 
on the properties of the stellar winds, which are in turn strongly 
linked to the total luminosity of the star, since these winds are 
radiation-driven. Knowing and understanding this Lx - ifiOL re- 
lati on is thus of prim e importance. 

iBerghofer et alj (fT997) investigated this scaling law for a 
large sample of objects (237 OB stars). Their study relied on 
data from the ROSAT All-Sky Survey (RASS), which provided a 
rather homogeneous dataset, and represented the first large-scale 
survey of massive stars in the X-ray range. The derived Lx-L^ol 
relation exhibits a large scatter (cr of 0.4 on a logarithmic scale, 
or a factor 2.5), and breaks down below Lboi ~ lO"*^ erg s"', i.e., 
for mid and late B stars. 

In recent years, new studies investigating this so-called 
'canonical' relation have been performed, using detailed 
observations of rich open clusters and associations, no- 
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tably NGC 6231 (X MM-Newton. ISaiia et al.1 l2006h . Carina 
OBI (XM M-Newton, lAntokhin et all 1200^ Westerlund 2 
(Chandra, iNaze et all l2008b . and Cyg OB2 (Chandra, 
lAlbacete Colombo et alj|200 7). For the clusters with the most 
well constrained stellar content (NGC 6231 and Carina OBI), 
the derived Lx - Lboi relations are far tighter (dispersions of only 
40%) than in the RASS analysis. 

The difference between the RASS and XMM-Newton results 
may be caused by two factors. First, the datasets and data han- 
dlings are not directly comparable: on the one hand, the RASS 
study is based on a large sample of stars with approximately 
known properties (count rates, spectral type); on the other hand, 
the recent studies have considered a small sample of well-known 
stars with precisely known properties (derived from X-ray spec- 
tra and optical monitoring). Second, the nature of the samples 
also differ: mixed field and cluster stars on the one hand, and 
a homogeneous stellar population from a single cluster on the 
other Metallicity and age could indeed affect the star's X-ray 
properties. 

To determine which factor is the most important, I decided 
to investigate the high-e nergy properties of a lar ge sample of hot 
stars, in a similar way as IBerghofer et al J d 1 9971) but by employ- 
ing the higher sensitivity and resolution of XMM-Newton. At 
the present time, two XMM-Newton surveys are available: the 
2XMMi catalog and the XMM slew survey (XMMSLl) catalog. 

In Sect. 2, 1 describe the main properties of the 2XMMi cat- 
alog, in Sect. 3 the sample of hot stars detected in the 2XMMi 
catalog, and in Sect. 4 the sample of hot stars detected in the 
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XMMSLl catalog. Sections 5 and 6 present the discussion and 
conclusion, respectively. 



2. The 2XMM\ catalog 

The XMM-Newton obse rvatory was launched in December 
1999 (iJansen et al 1 12001h . Its main instruments are three non- 
dispersive cameras, EPIC MOSl, MOS2, and pn ("Striid er et all 
[200 1; Turner et al., 2001), and two grating spectrographs, RGSl 
and 2 (Iden Herder et al.U200ll) . The EPIC cameras are sensitive 
to the 0. 15-12. keV band, whereas the RGS studies instead the 
0.35-2.5 keV band. While the RGS provides the highest spectral 
resolution {R-30Q at 1 keV), medium-resolution spectroscopy 
(R-IO at 1 keV) can also be obtained using the intrinsic energy 
resolution of the CCDs that make up the EPIC cameras. In addi- 
tion, time series can be derived from EPIC data, enabling vari- 
ability studies of the detected sources. 

Since pointed observations often detect many serendipitous 
sources in addition to the main target, the XMM Survey Science 
Centre (SSC) consortium compiled a catalog of serendipitous 
X-ray sources using 41 17 XMM-Newton archival datasets. This 
was developed to help fully exploit the capabilities of XMM- 
Newton, by means of a dedicated, homogeneous processing 
of all available data. The second main release of this cata- 
log, called the 2XMM catalog, was published on-line in 2007 
dWatson et al.L 12009') and an incremental version (2XMMi), 
more complete, was made available in the following year. A slim 
version of the 2XMMi catalog, containing the most important 
information, can be queried using VizierQ while the full version 
can be downloaded from the Vizier ftp site 0. Detailed informa- 
tion (images, spectra, time series) as well as a user guide are 
available at the SSC websiteB 

In this catalog, each source has a unique identifier based 
upon the lAU naming convention: it starts by 2XMM or 2XMMi, 
followed by the truncated position of the source at equinox 
2000.00 (JHHMMSS.S+DDMMSS). If the source was detected 
in several exposures, additional entries are available, each be- 
ing identified by a unique identifier (DETID) preceding the lAU 
name. More than 300 columns follow the name. They list the de- 
tails of the observations (observation identifier, filters used, ex- 
posure times,...), as well as the basic source properties (count 
rates, hardness ratios) in each instrument and for several en- 
ergy bands. Fluxes are also available but are not considered 
here because the count-rate-to-flux conversion used in the cat- 
alog considers a power-law spectrum typical of background 
AGNs, which indeed does not apply to hot stars. In this pa- 
per, the count rate in the 0.5-4.5 keV energy band (aka "band 
9") is preferentially considered since this rather broad band was 
used for a separate run of the SAS analysis tasks, which en- 
sured a better handling of the error values. Fortunately, it is also 
where hot stars emit most of their X-ray flux and the signal- 
to-noise ratio reaches its maximum. The four hardness ratios 
used in the 2XMMi catalog are defined as HR - where 
S=0.2-0.5keV and H=0.5-1.0keV for HRu S=0.5-1.0keV and 
H=1.0-2.0keV for HR2, S=1.0-2.0keV and H=2.0^.5keV 
for HR^, and S=2.0-4.5keV and H=4.5-12.0keV for HR^. 
Undetermined values are quoted as "NULL", which is the case 
for non-detections in one instrument (either because it was not 



' http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=IX/40 
^ ftp://cdsarc.u-strasbg.fr/pub/cats/IX/40 

^ http://amwdb.u-strasbg.fr/2xmmi/catexpert and 
|http://xmmssc-www.star.le.ac.iik/Catalogue/2XMMi/UserGuide-xmmcat 
respectively 



turned on or because the source fefl on a CCD gap). Time series 
and spectra are extracted for each detection of a given source if 
it presents at least 500 EPIC counts. 

It must be underlined that this catalog should not be consid- 
ered as a complete, homogeneous all-sky survey. First, the entire 
set of the regions observed by XMM-Newton cover only 420 
square degrees, or 1 % of the entire sky, and this coverage is very 
patchy. Second, these regions were not randomly chosen, but re- 
quested by different Pis for their specific needs, which means 
that the catalog is certainly biased: for example, about 65% of 
the observed regions have large Galactic latitudes \b\ > 20° be- 
cause of the strong interest in extragalactic fields. Finally, the 
exposure times range from shorter than 1000s to 130000s, and 
therefore the detection limit varies greatly from one area to an- 
other. However, as imperfect as it might be, it nonetheless con- 
stitutes the most sensitive X-ray catalog available at the present 
time in the observed regions, surpassing the RASS in every way. 
With -220000 unique detections, the 2XMMi catalog is the 
largest X-ray catalog ever produced: it contains twice as many 
discrete sources than previous surveys, notably the RASS. 



3. Finding hot stars in the 2XMM\ catalog 

To find the hot stars detected by XMM-Newton, it is necessary 
to correlate the 2XMMi catalog with a catalog of hot stars. To 
this aim, the Reed catalog of hot stars jReedl l2003l) was chosen 
because it is one of the most complete for this type of objects. 
The catalog version of January 2009 was kindly provided by its 
author. It contains >19 000 objects, mostly O and B stars, plus a 
few emission-line stars (WRs, T Tauri, HAe/Be,...). 

The correlation between these two catalogs was measured 
twice, once using the preferred "source name" in the Reed cat- 
alog and a second time using the coordinates provided by Reed 
(whenever Simbad was unable to resolve the provided name). 
The search radius was fixed to 5". This value represents a com- 
promise. On the one hand, the XMM-Newton PSF has a FWHM 
of ~5" on-axis but this greatly increases towards the edges of 
the field-of-view. The position uncertainty is expected to be a 
fraction of the PSF, and 96% of the detections indeed have a po- 
sitional error of <5". On the other hand, residual shifts can exist 
between X-ray and optical observations and the optical coordi- 
nates are certainly not 100% accurate. 

If two (or more) X-ray sources were detected within the cho- 
sen correlation radius, the closest one was kept only if the sec- 
ond closest source was more distant by at least 1" (to avoid spu- 
rious multiple detections). In addition, extended sources were 
discarded, since these objects can either be true diffuse X-ray 
sources, sub-areas of extended emission, or spurious, combined 
detections of several stars in a dense cluster High-mass X-ray bi- 
naries recorded in Simbad were also discarded, since this paper 
considers only the intrinsic emission of hot stars, and not their 
interaction with compact objects. Finally, a few detected sources 
were found to be listed twice in the Reed catalog: in these cases, 
only the smallest number used as an identifier in the Reed tables 
was listed in the correlation catalog 

The 2XMMi catalog provides a "quality flag" S , which is 
set to be zero if the X-ray source has no problem and to 1^ if 
it might be spurious. Experience showed that hot stars in clus- 
ters are generally detected with a non-zero quality flag, because 



This is the case of ALS9458=ALS4879, 9465=4887, 9468=4889, 
9477=4894, 9500=4912, 9506=4914, 9511=4918, 9526=4923, 
,Mi5fl6=4929, 9537=4930, 9604=5025, 9609=5039, 9621=5046, 
' 15858=1855, 21108=11417. 
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of the presence of numerous neighbours. To be as complete as 
possible, no selection was made based on the quality flag and 
all detections were kept. However, the higher quality objects are 
separated from the others in the discussion and tables provided 
below. To check whether spurious assoc iations were a problem, 
the approach of iBerghofer et al.l (Il996h was used. The number 
of background objects is expected to scale with the number of 
detected hot stars, the area of detection (5" radius), and the X- 
ray source density (about 500 src per square degree, compared 
to the 1.5 src per square degree of the RASS): for our data, only 
one spurious source is expected. Inspecting all bright sources 
(i.e., those with available light curves and spectra) identified 
only one really spurious detection, associated with ALS 4440 
(=HD3 16464), which is part of the PSF wing of the bright X- 
ray binary XTE J1751-305. 

The result of the correlation between the Reed and 2XMMi 
catalogs is presented in Table 1 (available electronically at the 
CDS). The first column provides the data used to query the cat- 
alog, i.e., either the prefered name in the Reed catalog or the 
Reed coordinates. It is preceded by a letter in case two stars 
correspond to the same X-ray source. The second column lists 
the Reed number of the star (ALS # in Simbad), while the third 
one provides the usual name of the source (HD, BD, CPD,... 
identifiers). The fourth column indicates the spectral type of the 
source, which was chosen to be the most recent from the Reed 
catalog, except for the Wolf-Rayet s tars (WRs) where the clas- 
sification of the Vllth WR catalog (Ivan der HuchA 1200 ll) was 
used. If no type is available in the Reed catalog, the spectral 
type from Simbad is reproduced here, preceded by a ~ to eas- 
ily distinguish these sources. The fifth column corresponds to 
a flag set to be Y if the source was detected as a binary. A 
star was classified as binary if either the spectral type avail- 
able in the Reed catalog indicates the presence of a compan- 
ion, or if the star i s a kn own bina ry in the 9th Binary cata- 
log (Pourbai x et"an |2004) and/or in iOiesI ('20b3). The next few 
columns indicate the V magnitude, B-V colour, interstellar ab- 
sorption column, and bolometric flux. The photometric data cor- 
respond to the most recent values in the Reed catalog or, if un- 
available, to the Simbad values (again flagged with a ~). For 
the WRs, the interstellar a bsorption column and the bolomet- 
ric fluxes were taken from Oskinoval (l2005h for WN stars. For 
O and B stars, these characteristics were calculated whenever 
a precise spectral type (i.e., not only 'O star' or 'B star') and 
magnitude values were known. The interstellar column was de- 
rived from B ohlin's formula (A^^ = 5.8 x 10^' x E(B - V) cm"^, 
iBohlin et aljri978i) . where the colour excesses were calculated 
from the difference between the observed colours and the in- 
strinsic ones. Because of the scarcity of accurate distances, bolo- 
metric fluxes were preferred to bolometric luminosities. They 
were calculated using the usual formulae, yielding log(/BOL) - 
-4.61 - where /bol is in ergcm"^ s"'. Intrinsic 

colours and bolometric corre ctions (BC) were taken , for the con- 
sid ered spectral types, fro m iMartins & Pled (l2006h for O stars 
and lSchmidt-Kalei](II982h for B stars. In the case of binaries, the 
intrinsic properties (colour, BC) of the primary star were used. 
The subsequent columns correspond to the X-ray source prop- 
erties, i.e., the distance between the X-ray source and its optical 
counterpart, the official lAU name of that counterpart, list of the 
filters used for each EPIC instrument (pn, MOSl, and MOS2 
in this order), the total count rates and their errors in the 0.5- 
4.5 keV band for each instrument (with pattern 0^ for pn and 
0-12 for MOS), the four hardness ratios for each instrument, the 
total detection likelihood, the overall quality flag (see previous 
paragraph), two variability flags (see Sect. 3.1.2), the number 



Table 2. Classification of the optical sources detected in the 
2XMMi catalog. The number in parentheses yields the number 
of sources with available EPIC spectra. 



Quality flag 


O stars 


B stars 


WR stars 





35 (16) 


91 (22) 


7(6) 


i=0 


98 (62) 


70 (28) 


9(8) 



Table 3. Percentage of known binaries amongst the detected hot 
stars. 



Quality flag 


stars 


B stars 


WR stars 





17% (6/35) 


5% (5/91) 


57% (4/7) 


^0 


22% (22/98) 


4% (3/70) 


89% (8/9) 



of individual detections, the unabsorbed X-ray flux in the 0.5- 
10. keV band (see Sect. 3.2), and finally the log Lx/Lbol value. 

If several observations were available for a given source, a 
mean count rate was calculated for each instrument-i-filter con- 
figuration. EiTor bars were compared for different filter choices 
of a given instrument; the instrumental configuration corre- 
sponding to the smallest error bar was considered as the most 
suitable one and reproduced in Table 1 . 

Note for individual objects: ALS 4440 is spurious (see 
above); the counterpart of ALS 4592 (=CPD -24°6140) is real 
but affected by straylight, so that its properties might not be to- 
tally uncontaminated; ALS 20163 (=West 1 30) shares its X-ray 
counterpart with a sgBe star. 

3. 1 . Basic properties of tiie detected objects 

From the coiTelation performed above, 310 stars (with a known 
O, B, or WR spectral type) were found to have an XMM-Newton 
counterpart within 5". Of these, 133 have a zero quality flag and 
correspond to 132 unique 2XMMi sources (i.e., two sources are 
found within the PSF of one X-ray counterpart). The remain- 
ing 177 objects are considered as potentially spurious within 
the 2XMMi catalog; they coiTespond to 173 individual 2XMMi 
sources, i.e., there are four pairs of stars sharing the same X-ray 
counterpart. 

The biases inherent to the generation of this 2XMMi cata- 
log have already been pointed out in Sect. 2 and the conclusions 
found from the performed correlation should indeed be inter- 
preted with caution. Nevertheless, I examined the basic proper- 
ties of the detected objects. While they may not be totally rep- 
resentative of the entire population of hot stars, these detections 
certainly represent the largest sample of such stars studied by 
modern X-ray observatories. Note that the remainder of the pa- 
per focuses mainly on the detailed properties of OB stars, the 
spectral analysis of WR stars being defeiTed to a future paper. 

3.LL Distribution 

The repartition of the detected stars amongst the main spectral 
types is given in Table |2] and shown in Fig.[T] For sources with 
a good quality flag, the majority display a B spectral type, as 
one could have expected based on the incidence of these lower- 
mass objects compared to O and WR stars. When potentially 
spurious sources are included, the situation appears to change 
because the O stars now dominate. This is probably caused by 
the combination of two facts: (1) the O stars are brighter in X- 
rays, and (2) these "spurious" sources are found mostly in clus- 
ters, where only bright O-type stars are easily detected against 
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Fig. 1. Distribution of the spectral types and luminosity classes of the O and B stars detected in X-rays (when detailed classification 
is available; for binaries, only the classification of the primary was considered). The middle panels correspond to X-ray sources with 
a good quality flag (5'=0), while the upper panels show the results for "potentially spurious" sources (5'!=0). The corresponding 
distributions for the entire Reed catalog are shown in the lowest panels, where the dotted red line indicates the scaled distribution of 
all detected objects (zero and non-zero quality flags, with Poissonian errors). Note that a luminosity class of 6 indicates stars with 
known spectral type but undefined luminosity class. 



the overall background produced by the combined X-ray emis- 
sion from PMS stars, background objects, and diffuse emission 
within the cluster. 

Interestingly, it must be noted that, excluding WR systems, 
binarity does not seem to play an important role in the probabil- 
ity of X-ray detection (see Table [3] only ~20% of the detected 
O stars and ~5% of the detected B stars are known binaries). 
This result might first seem at odds with previous conclusions 
based on older data (e.g., Einstein observations in ChlebowskH 
119891) . Although one cannot exclude the possibility that part of 
the binary population remains undiscovered (especially for the 
B stars), this bias cannot explain the new result since the older 
data were affected by similar uncertainties; besides, it would in- 
deed be unexpected that XMM-Newton targeted a specific area 
of the sky with a low number of binaries. More importantly, 
it must be recalled that the preferred detection of binaries was 
then explained by the presence of an additional X-ray emission 
produced by wind-wind collisions. However, it was shown that, 
in sensitive observations of an entire hot-star population, only 
a small fraction of the massive binaries display strong wind- 
wind collisions capable of emitting X-rays: only a f ew systems 
are th us truly overluminous in the 0.5-lOkeV range (ISana et al.L 
|2006|) . This is linked to the X-ray excess, when detected, be- 
ing possibly dependent on the bolometric luminosity and the 
orbital parameters (Chlebowski & Garmanv, 1991; Linde r et"an 
I2OO6 ). Therefore, there may be two different causes of the appar- 
ent discrepancy between more recent observations and the first 
Einstein results: the sample size (a few bright, peculiar objects 
compared to a larger sample of many different cases), as well as 
improvements in the knowledge of the physical parameters (e.g., 
with respect to stellar multiplicity, especially the detection of 
faint companions) and in the X-ray instrumentation (higher spa- 



tial/spectral resolution and higher sensitivity, including at lower 
energies than Einstein, i.e., where the wind-wind emission is 
less prominent than at higher energies). 

Finally, the distribution of the spectral types and luminos- 
ity classes is shown in Fig. [T] Amongst O-type stars, the more 
numerous 06-9.7 stars constitute most of our sample. This is 
also true for main-sequence systems, which are more numerous 
in general, this study being no exception. The O-star distribution 
well reflects that of the initial stellar catalog, as can be seen in the 
bottom panels of Fig. [T] However, there is a lack of detections 
amongst late B-type stars, while earliest B-stars as well as giant 
B-stars appear preferentially detected. At first sight, the former 
results appear in phase with the current wind-shock model: one 
does expect the X-ray luminosity to drop (hence the detection 
rate) as the stellar wind weakens and then disappears. 



3.1.2. Short- and long-term variations 

Table |4] summarizes the results of the variability study. The 
short-term variability was directly analyzed during the 2XMMi 
processing. The first variability flag quoted in Table 1 corre- 
sponds directly to the variability flag of the 2XMMi catalog, i.e., 
it is defined to be 1 if the source was found to be variable with 
a significance level of 0.001% following a.x^ test performed on 
the time series of the individual exposures. For hot stars, as could 
be expected, very few objects vary within one exposure (<10%). 
Three sources display a flare, which is quite typical of low-mass, 
pre-main-sequence stars: HD37016 (B2.5V), HD33904 (B9III), 
and HD37479 (crOriE, B1/2V, whose flare decay has a rather 
long time constant). In addition, HD120991 (B2Ve) exhibits a 
clear increase of its count rate during the observation. For the 
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Table 4. Number and percentage of variable X-ray counterparts. When comparing several exposures, the percentages noted for 
non-tested objects were calculated with respect to the entire sample, whereas those of the variable sources refer only to the tested 
subsample. 



Q. flag 


Sp. Type 


Within one exposure 


Between exposures 








total 


single 


binaries 


total 


single 


binaries 


not tested 

















4 (25%) 


3 (25%) 


1 (25%) 


19 (54%) 





B 


4 (4%) 


4 (5%) 





7 (32%) 


7 (33%) 





69 (76%) 





WR 











2 (50%) 


1 (50%) 


1(50%) 


3 (43%) 




O 


4 (4%) 


3 (4%) 


1 (5%) 


44 (69%) 


32 (67%) 


12 (75%) 


34 (35%) 


to 


B 


4 (6%) 


4 (6%) 





14 (67%) 


12 (63%) 


2 (100%) 


70 (49%) 




WR 











6 (86%) 





6 (100%) 


2 (22%) 



other stars, there is no obvious flare, and the variations cannot 
easily be differentiated from those of the background signal. 

Long-term variability, for example between exposures, can 
also occur, but the 2XMMi catalog does not check for its pres- 
ence. A second variability flag was thus calculated for Table 1: 
it is set to be 2 if there are not enough exposures (0 or 1 observa- 
tion for all combinations of filters/instruments), to 1 if a test 
detected variations in the count rate with a significance level of 
1%, and otherwise to 0. As is obvious from TablelH this type of 
variability is far more common in hot stars, with few differences 
between single and binary OB stars. Except for the putative pres- 
ence of colliding-wind binaries or magnetic wind confinement, 
the cause of these variations remains unknown. 

3.2. Spectral fitting for tlie X-ray brigiit objects 

About half of the detected hot stars (59% of the O stars, 31% of 
the B objects, and 88% of the WRs) have enough counts within 
a single exposure to have their spectra automatically extracted 
by the 2XMMi processing (Table |2]l. Many of these were ob- 
served several times: the individual spectra were merged, instru- 
ment by instrument, only if the source was non-variable. In total, 
332 spectra of OB stars were finally fitted. They were individu- 
ally fitted, in the 0.3-10. keV energy range, within Xspec vll.2.0 
with absorbed multi-temperature thermal plasma models of the 
type wabs{IS M)xwabsx('^ mekaU) with i equal to 3 at most and 
solar abundances assumed. The first absorbing column was fixed 
to the ISM column derived above (see Sect. 2) if known or to 
otherwise. The second absorbing column allows for additional, 
local absorption. A word of caution should be added: the soft 
and hard thermal components might not be formed at the same 
location inside the wind and might thus be affected by different 
local absorptions. However, we restrict ourselves to one common 
local absorbing column because, for the majority of the stars in- 
vestigated in this paper, the quality of the data does not justify 
the use of a model of the type wabs{IS M) x (2 wabs, x mekaU). 
The number of optically-thin plasma components is the mini- 
mum number of single-temperature components needed to pro- 
vide a good fit (x^ ~ 1 or < 2) to the considered data. Generally, 
the IT-models are used only when there are too few counts, i.e., 
only the main component is distinguishable. When the signal- 
to-noise ratio is reasonably high, the sum of two thermal com- 
ponents is requested to provide a good fit. 

Table 5 (available electronically at the CDS) provides the re- 
sults of these fits. The first column indicates the Reed number of 
the objects, the second column their usual name, and the third 
column the interstellar absorbing column (N^^, in cm"^). The 
fourth and fifth columns provide the value of x^ and the number 
of degrees of freedom. The next four columns list the absorbed 
fluxes (inergs-i cm^^) in the 0.5-10. keV, 0.5-1. keV, 1.-2.5 keV, 



and 2.5-10. keV energy bands. They are followed by 4 similar 
columns providing the unabsorbed fluxes (i.e., the X-ray fluxes 
dereddened by the insterstellar component, if known, or other- 
wise by the fitted absorption column). The additional absorption 
is given in Col. 14 (A^h, in 10^^ cm"^), followed in the next two 
columns by the lower and upper limits of the 90% confidence 
interval. The temperature (kT in keV) and normalization factors 
(norm in cm"^) of each mekal component in turn (in order of 
increasing temperatures) are then given, together with their 90% 
confidence interval. Finally, the last column provides additional 
remarks, such as the observation identifiers (if several exposures 
were available) or, in a few cases, any modification of the gen- 
eral fitting scheme (e.g., non-solar abundances). Undetermined 
values are again quoted as "NULL". Note that a comparison was 
made between our results and published analyses of some hot 
stars' spectra ( Of?p stars. Naze et al. 2004, 2007, 2008; flCar, 
Naze & Rauw i 120081: N GC6231. iSana et al. 2006; HD1 68112 
De Becker et a l."2004'; Cyg OB2 #8A, ,De Becker eTall 120061; 
Carina OBI, lAntokhin etai]|2008|) and a good agreement was 
found. It must also be emphasized that these fits provide an ade- 
quate idea of the flux distribution for the spectral resolution and 
signal-to-noise ratio considered, but they should not be taken at 
face value (especially in the case of non-solar abundances), since 
the true properties of the hot plasma can differ (e.g., continuous 
distribution of thermal components versus a few discrete tem- 
peratures). 

Figure |2] graphically illustrates the results. In terms of the 
absorption, the interstellar components appear larger for O stars 
than for B stars, as expected (O stars being brighter and hence 
detectable from further away). In addition, a significant, addi- 
tional absorption is needed when fitting O-star spectra. Since 
its distribution does not correlate with the interstellar absorp- 
tion histogram, this effect should be considered to be signifi- 
cant. Concerning temperatures, different cases corresponding to 
different number of mekal components should be envisaged. If 
only one thermal component was fitted, its temperature is usu- 
ally quite low, 0.2 or 0.6 keV, for O stars, while the distribution 
is much flatter for B stars, temperatures in half of the cases being 
below 1 keV and in the other half above. If two thermal compo- 
nents are fitted, the first temperature is generally quite soft, with 
a frequent ambiguity between 0.2 and 0. 6 keV ( for previous re- 
ports, see e.g., the case of HD148937 in lNaze et al. 2008). The 
second temperature is about 2 keV for O stars and B stars, but 
has a far larger dispersion than in the latter case. If three ther- 
mal components were fitted, the first temperature was always 
very low (often at the lowest limit of the permitted range within 
Xspec, 0.08 keV), the second temperature again is 0.2 or 0.6 keV, 
while the third one is generally 2 ke V. The normalization factors 
of the additional components usually have lower values than that 
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Fig. 2. Distribution of the spectral properties (temperatures kr„ absorptions and A^h, ratio of normalization factors 

norm,/normi) for the whole dataset, the O-type stars only and the B-type objects only. Note that, when several fits were made, 
an average value was used in the computation of the histogram and that bins at the rightmost limits of the plots correspond to all 
values larger than this limit. 



of the first component, except in the case of 2-temperature fits to 
B stars, where the harder component clearly dominates. 

3.2.1. X-ray fluxes 

Spectral modelling directly provides the X-ray fluxes of each 
fitted object (see Table 5), which are simply reproduced in Table 
1. In the case of a varying source observed several times, the 
X-ray flux reported in Table 1 is the average flux. 

However, the objects in the fainter half of the sample do not 
have spectra, and any estimate of their flux must rely only on 
their count rates. To convert count rates into fluxes, the aver- 
age/the most probable 2T-model was used. For O-type stars, it 
has a local absorbing column of ~4xl0^'cm"- and tempera- 
tures of 0.2 and 2keV, which is represented well by the spec- 
trum of HD168112 (ALS 4912, second observation): the latter 
was therefore used as a reference. For B-type stars, the aver- 
age properties are a zero local absorbing column and tempera- 
tures of 0.2 or 0.6 and 2 keV, and a good reference in this case 
is HD37040 (ALS 14653). The spectral properties of these two 
reference stars (A^H,kri 2, flux ratio at 1 keV of 0.158 and 1.46, 
respectively) were used as input for PIMM^ The unabsorbed 
fluxes (i.e., dereddened only for the interstellar component) were 
calculated for a unit count rate, one of the two cameras, one of 
the three filters, and a range of interstellar absorbing columns 
(from to 2.75 xlO^' cm"^). The resulting conversion factor had 
to be corrected by a factor of 0.70 because PIMMs considers 
count rates in a region of 1 5" radius whi l e the 2XMMi catalog 
provides the full count rate (IWatson et al.L |2009I FI For each star, 
the X-ray flux was calculated by taking into account its own in- 
terstellar absorption (see Sect. 2), no value was thus quoted if 
this absorption is unknown. The results were checked against 
the bright stars and good agreement was found - in 80% of the 
cases, the flux derived from the count rates is within a factor of 2 



' http://heasarc.gsfc.nasa.gov/Tools/w3pimms_pro.html 
* In addition, PIMMs considers only counts with a zero PATTERN 
keyword, while the 2XMMi counts are for PATTERN keywords rang- 
ing from to 12 (MOS) or to 4 (pn). However, as most of the recorded 
counts have a zero PATTERN, this uncertainty is considered to be neg- 
ligible compared to other sources of error. 



of the flux derived from the spectral fits. The derived unabsorbed 
X-ray fluxes for the faintest stars can be found in Table 1 . 

Once the fluxes are known, the log/x//BOL ratios, equiva- 
lent to the logLx/i'BOL ratio, were derived. They are shown in 
the last column of Table 1 . Note that, for the five pairs of stars 
sharing the same 2XMMi counterpart, the ratio is calculated for 
the total bolometric flux. The unweighted average log/x//BOL 
ratios (for all stars but the pairs) were then calculated, together 
with their dispersion. The results are listed in Tableland shown 
in Fig. 1^ 

4. Detections of hot stars in tlie XMM-Neivfon Slew 
Survey 

Between pointed observations, XMM-Newton slews between 
different sky areas. During many of these slews, the telescopes 
remain open and thus perform a limited survey. The analysis of 
the EPIC-p n data from 218 slew s detected more than 2600 X- 
ray sources dSaxton et 311120081) . They constitute the XMMSLl 
catalogue, which can be queried using Viziefl 

As for the 2XMMi catalog, the results of the slew survey 
were cross-correlated with the preferred names and coordinates 
extracted from the Reed catalog, again for a correlation radius of 
5", yielding 13 detections. Extended X-ray sources and sources 
already belonging to the 2XMMi catalog were discarded from 
the list, as well as wrong identifications (i.e., a hot star in the 
vicinity of an X-ray source, but where the hot star is not the 
"official" counterpart listed in the slew survey catalog). In addi- 
tion, counterparts were automatically identified for slightly more 
than half of the X-ray sources of the XMMSLl : their nature was 



' It might at first appear surprising that the detected O and B stars dis- 
play similar bolometric fluxes. It must however be kept in mind that the 
log Lx/LgoL ratios for these stars are rather similar: if the X-ray sources 
constitute a flux-limited sample, the bolometric fluxes of the detected 
objects will also appear quite similar. Indeed, this conclusion does not 
apply to bolometric luminosities: as shown in Fig. |2] the O stars of 
our sample are more absorbed (hence intrinsically much brighter after 
correcting for the ISM absorption), and much more distant, than the B 
stars. 

** http://vizier.u-strasbg.fr/viz-bin/VizieR?-source=I/A+ A/480/6 1 1 
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Table 6. Average log f^""''y /bol ratios, with their dispersions. 



O stars 



B stars 



From fits only 
0.5-10 keV 

0. 5-1. keV 

1. -2.5 keV 
2.5-10 keV 



-6.45 ±0.51 

-6.69 ± 0.45 
-7.02 ± 0.58 
-8.03 ± 1.32 



-6.27 ± 0.98 

-6.73 ± 0.69 
-6.82 ± 1.07 
-7.48 ± 1.89 



From fits and conversion of the count rates (0.5-10 keV band) 



all 


-6.66 


±0.50 


-6.59 = 


tO.83 


5=0 


-6.80 


±0.35 


-6.81 = 


tO.74 


5 to 


-6.61 


±0.54 


-6.34 = 


tO.87 


binaries 


-6.47 


±0.52 


-7.34 = 


tO.62 


singles 


-6.72 


±0.49 


-6.55 = 


tO.82 


supergiants 


-6.75 


±0.36 


-6.82 = 


tO.88 


giants 


-6.74 


±0.33 


-7.13 = 


tO.85 


main sequence 


-6.67 


±0.49 


-6.33 = 


tO.77 


Carina 


-6.54 


±0.56 






NGC6231 


-6.71 


±0.22 






NGC6604 


-6.48 


±0.48 






CygOB2 


-6.97 


±0.40 







checked using Simbad, which inferred that they were 34 X-ray 
emitters associated with hot stars, sometimes more distant than 
5". Of these, those stars not belonging to the Reed catalog or 
those already found in the 2XMMi catalog were discarded. This 
two-step trimming process finally resulted in the detection of 16 
additional X-ray sources, presented in Table [T] where the total 
count rates are given for the pn instrument with the medium fil- 
ter and in the 0.2-12. keV band, and the DL column refers to the 
individual detection likelihoods. 

As for the 2XMMi data, the optical properties of the star 
(e.g., interstellar reddening, bolometric flux) were derived from 
their known characteristics (e.g., spectral type, magnitudes). The 
count rates were then converted into unabsorbed fluxes in the 
0.5-10. keV band, in a similar way to the method described in 
Sect. 3.2. 1, to derive the Lx - Lbol ratio. 

5. Summary and discussion: the X-ray properties of 
OB stars 

Since the slew survey provided only a few additional sources 
with poorly constrained X-ray properties, the following conclu- 
sions rely mostly on the detailed analysis based on the 2XMMi 
catalog. 

Since the XMM-Newton observations cover only 1 % of the 
sky, a first question arise about the possible detection biases. 
For the detected O stars, the distribution of their spectral types 
and luminosity cl asses i s fully comparable to that of the full hot 
star catalog from Reedl (l2003i) . so that no bias whatsoever is ex- 
pected. This is not the case for the detected B stars, where a lack 
of late-type objects, as well as an excess of early-type stars and 
giant objects, are found. Although the consequence of selecting 
specific celestial regions should not be disregarded, it is proba- 
ble that this effect is real: since the very first X-ray detection of 
hot stars, it has been shown several times that early-B stars are 
brighter, hence more easily detected, than late-type objects. 

Variability studies demonstrate that short-term changes are 
quite rare. When they exist, they generally appear flare-like, as 
found for PMS stars. It is as yet unclear whether the emission 
originates in a PMS object (true companions or line-of-sight co- 



incidences), or if a magnetic phenomenon similar to that oper- 
ating in PMS stars is at work in those objects. Stochastic varia- 
tions such as those expected from collisions between wind-shells 
or large clumps are not seen. Long-term variations appear to be 
far more common in hot stars since they affect about half of the 
cases. They could well be caused by e.g., wind-wind collisions in 
unknown binaries, but will require an additional, in-depth anal- 
ysis before a definite cause can be identified. 

The spectral fits have enlighted the differences between O 
and B stars. The X-ray spectrum of O-type objects requires an 
absorption in addition to the interstellar component. This is not 
an effect of the uncertainties in the latter, which is quite well 
constrained since the intrinsic colour of O stars is quite simi- 
lar regardless of subtype or luminosity class (altho ugh it is true 
that th e Nh - E(B - V) plots display some scatter: iBohhn et akl 
Il978h . It thus appears intrinsic to these stars, most probably re- 
lated to the presence of their stellar winds. Indeed, B-type stars, 
whose winds are much weaker, need not have any additional ab- 
sorption. Another difference lies in the temperatures needed to 
fit the plasma emission. On average, O-type stars exhibit a soft 
emission, well fitted by a thermal component with a temperature 
of 0.2 or 0.6 keV, with possibly a faint hard component (typi- 
cally at 2keV). On the other hand, the spectra of B-type stars 
are harder, well fitted by hot thermal components (either a sin- 
gle one at about 1 keV, or the sum of a faint warm one with a 
temperature of 0.2-0.6 keV and a brighter hot component with a 
temperature of 2keV). These hot components are not expected 
in the usual wind-shock model, and require the presence of either 
a PMS star along the same line-of-sight (physical companion or 
not) or of exotic magnetic phenomena. Analyses of multiwave- 
length follow-up observations would be required to settle this 
question. 

Finally, the Lx - L^ol ratio was also investigated. 
Unsurprisingly, the dispersion is much lower for O stars than 
for B stars, and for the soft or medium energy bands than for 
the hard band (see e.g., similar result s for specific clusters in 
ISana et all 120061: lAntokhin et ani2008l) . The relation for the B 
stars also appears shallower and might have a more compli- 
cated parametrization (e.g., log(Lx) -ax log(LBOL) + b rather 
than log(Lx/LBOL) = b) but this cannot be tested here because 
of a lack of accurate distances. For O-type stars, the disper- 
sion in the Lx - L^ol relation is about 0.35-0.5 dex, or a fac- 
tor of 2-3, a value similar to that found in the RASS analysis 
(Berghoferet al., 1997). Note that this dispersion is not mostly 
caused by uncertainties in the fluxes, since the fluxes are derived 
from precise fits, and those calculated from a general conversion 
of the count rates are similar in value. It thus seems that the dis- 
persion measured by Berghofer et al. is real, and the much tighter 
correlation found for some clusters (e.g., NGC6231, Sana et akl 
2006 and Table ^ needs to be explained by reasons other than 
a different treatment of the data. Although the dispersions are 
somewhat large and the samples quite small, the Lx -Lbol ratios 
found for different clusters indeed appear to differ marginally. 

It must also be underlined that no significant difference was 
found between the Lx - Lbol ratios of objects with various 
quality flags, binary status, or luminosity classes, although we 
are only just below the Icr limit in the case of B-type stars. 
As discussed above (Sect. 3.1.1), binaiity now appears less im- 
portant than advocated in the past. One possible explanation 
for the weak impact of wind-wind interactions could be radia- 
tive inhibition, where the radiation of one component deceler- 
ates the wind o r iginat ing in its companion (Gayley et al., 1993: 
lAntokhin et all |2004 . It should have more impact on short- 
period systems than on long-period binaries. Figure |4] shows 



Table 7, Additional X-ray detections of hot stars using the XMM slew survey. 



Reed name 


ALS# 


Usual name 


Sp. Type 


Bin.? 


V 


B-V 


d 

(") 


XMMSLl name 


DL 


Count rate 
ct s"' 


log(/BOL) 


funabs 

-2 -1 

erg cm s 


log(/x//BOL) 


HR3129 


14910 


HD 65818 


B1V+B3: 


Y 


4.51 


-0.18 


8.2 


J075815.2-491444 


9.7 


0.55+0.31 


-5.23 


l.Ole-12 


-6.77 


LSni+44 43 


11848 


BD +43 3913 


B1.5V 




8.91 


0.54 


6.4 


J212501. 9+442706 


15.8 


0.57±0.20 


-6.18 


3.68e-12 


-5.25 


LSE 99 


18922 


CD-35 9665 


0+ 




12.70 




4.3 


J143740.0-36132 




1.85±0.58 








LS 4635 


4635 


HD 165052 


06V((f)) 


Y 


6.87 


0.09 


3.2 


J180510.3-242354 






-5.47 






LS III +43 19 


11807 


HD 203064 


07.5111 


Y 


5.04 


-0.06 


5.0 


J2 11827.7+435646 


10.8 


0.65+0.25 


-5.06 


1.68e-12 


-6.72 


HR 3860 


16348 


HD 83979 


B5V 




5.07 


-0.14 


1.0 


J093353. 1-805629 




0.70+0.38 


-6.01 


l.lle-12 


-5.94 


HD 157832 


17479 


HD 157832 


B3Ve 




6.66 


0.00 


7.7 


J172754.0-470132 


20.1 


1.07+0.32 


-6.25 


2.57e-12 


-5.34 


LS V +35 1 


7841 


HD 24912 


07.5111(f) 




4.04 


0.00 


7.0 


J035857.3+354731 


11.4 


1.55+0.69 


-4.59 


4.58e-12 


-6.75 


HD 212571 


14728 


HD 212571 


Bllll/IVe 


Y 


4.80 


-0.18 


16.4 


J222517.0+012223 


18.6 


2.40+0.73 


-5.46 


4.41e-12 


-5.90 


HR6510 


15050 


HD 158427 


B3V 




2.83 


-0.11 


5.8 


J173149.8-495233 


15.4 


0.84+0.31 


-4.85 


1.59e-12 


-6.95 


HD 133242 


17476 


HD 133242 


B5V 




4.59 


-0.14 


3.5 


J150507.4-470303 




1.04+0.40 


-5.82 


1.66e-12 


-5.96 


Hbg 1339 


14476 


HD 155806 


07.5ine 




5.62 


-0.04 


8.6 


J171519.4-333246 


17.2 


0.72+0.25 


-5.27 


1.95e-12 


-6.44 


LS 4148 


4148 


WR93 


WC7+07/9 


Y 






15.7 


J172507.6-341110 


14.5 


0.53+0.19 








MCW 441 


14798 


HD46328 


BlIII 




4.33 


-0.24 


6.0 


J063 150.9-232509 


16.6 


1.40+0.42 


-5.34 


2.16e-12 


-6.33 


HR 3055 


14901 


HD63922 


BOIII 




4.10 


-0.20 


4.7 


J074913.9-462227 




0.96+0.27 


-4.98 


1.82e-12 


-6.76 


LS 4142 


4142 


HD319718 


03If* 




10.43 


1.45 


4.6 


J172443.5-341201 




0.69+0.22 


-5.11 


8.10e-12 


-5.98 
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Fig. 4. For the O+OB binaries with spectral information and de- 
fined period, log fx/ /bol ratio from Table |6] 



the variation in the log[Lx/LBOL] ratio with the binary period, 
when the latter is known precisely (this is the case for 23 sys- 
tems out of the 27 Oh-OB binaries having spectral information). 
For short-period systems (i.e., P < 35 d), there are hints of an 
overall shallow trend of increasing log /x //bol towards longer 
periods. Without a detailed modelling of each individual sys- 
tem, it is however difficult to assess whether the main cause 
of this trend is the radiative inhibition phenome non or the ex- 
pecte d Lx oc Mv^ variation for radiative systems (IStevens et al.L 
I1992I) . It has to be emphasized that the scatter in the points of the 
log /x //bol -period diagram is much larger than the mild "trend" 
mentioned above. For long-period systems (i.e., P > 35 d), the 
average log /x //bol appears smaller than for short-period bina- 
ries, which might reflect the well-known Lx°^ D^^ variation for 
adiabatic systems (Stevens et al., 1992), but the lower limit to 
the log /x //bol ratio appears to increase towards longer periods, 
which might be related to the smaller impact of the radiative in- 



hibition. However, the small number of systems in this range of 
periods (only 2!) prevents me from drawing solid conclusions. 

6. Conclusions 

This paper presents the results of the first global survey of hot 
stars with the highly sensitive X-ray observatory XMM-Newton. 
It relies on the 2XMMi and XMMSLl catalogs. About 330 stars 
were detected, re presenting a sam ple comparable in size to the 
RASS detections ("Berghofe r et al.L [19 96) but only covering 1% 
of the entire sky (to be compared with the full sky investigated 
in the course of the RASS). 

The derived properties of the hot stars confirm the results of 
the first preliminary data of the RASS and the recent in-depth 
investigations of a few clusters. The O stars have relatively soft 
spectra and exhibit a rather tight Lx - ^bol relation, although 
the dispersion is closer to that observed in the RASS survey than 
in a dedicated analysis of specific clusters. The B stars appear 
to have far harder spectra, with brighter emission for the ear- 
liest subtypes and a Lx - ^-bol relation with far larger scatter 
However, it must be noted that on average binaries do not ap- 
pear significantly brighter than single objects. 

The detection lists and the results of more than 300 spectral 
fits are made available to the community. 
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Fig. 3. Unabsorbed X-ray flux as a function of the bolometric flux. Lines indicate the average log /x//bol ratio from Tabled In the 
case of binaries, the luminosity class of the primary was used as a reference in the middle panels. 



